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Processes on different length scales affect the dynamics of chain molecules. In this work, we focus 
on structures on the scale of a monomer and investigate polyolefins, i.e. hydrocarbon chains with 
different small scale architectures. We present an exact enumeration scheme for the simulation of 
interactions and relative motion of two short chain sections on a lattice and employ it to deduce 
the probability for segmental motion for polymers of four different architectures in the melt. The 
probability for segmental motion is inversely proportional to the monomeric friction coefRcient and 
hence the viscosity of a polymer. Combining our simulation results with an equation of state 
for the thermodynamic properties of the polymers, we are able to make predictions about the 
variation of the friction coefficient with temperature, pressure, and small scale architecture. To 
compare our results with experimental data, we have determined monomeric friction coefficients 
from experimental viscosity data for the four polyolefins considered in this work. For temperatures 
well above the glass transition temperature, we find that our simple approach gives a good qualitative 
representation of the variation of the friction coefficient with chain architecture, temperature and 
pressure. 



I. INTRODUCTION 



The dynamics of chain molecules are affected by local interactions between individual chain segments as well as 
processes on the length scale of the whole chain and collective motions of chain segments. In this work, we investigate 
the effect of small-scale chain architecture on the dynamic properties of polymers. The polyolefins depicted in Fig. |l| 
are a good example for this eiiect. These polyolefins, all hydrocarbons with sum formula C„Il2n, differ considerably 
in their viscoelastic propertieglJTJ despite their chemical similarity. 

An important ingredient in theories for polymer dynamics is a friction coefficient ^ which is employed in coarse- 
grained models to describe small-scale effects on the dynamics of the system (cf. Ref. ^l]). In the Rouse model, for 
example, chain segments consisting of many monomers are represented by a single bead and spring. The Rouse 
friction coefficient Cr describes the-damping of the bead motion by the surrounding medium. The Rouse viscosity ryp; 
is proportional to and given bylij 

Na Rn 

where Na is Avogadro's number, p is the mass density, AI is the molecular mass, Rq is the mean-squared end-to-end 
distance of the chain, and Nji is the number of Rouse segments per chain. The Rouse model describes the viscosity 
of polymer melts of sufficiently short chains, i.e. for polymers with a molecular mass well iielow the entanglement 
mass Mq. For larger molecular masses, entanglement effects have to be taken into account.Ej While p and i?g vary 
with temperature at ctpnstant pressure, almost all of the temperature dependence of the viscosity is contained in 
the friction coefficient B Since the viscosity increases strongly as the glass transition temperature is approached, the 
friction coefficient C reflects collective motion as well as small-scale interactions. |— , 

The friction coefficient is inversely proportional to the probability for segmental motionJ13 a relationship that we 
are going to exploit in this work. In an exact enumeration procedure we perform lattice simulations of relative motion 
and interactions of two short chain segments, where the surrounding medium is represented in an average way. The 
collected statistics are evaluated as described in Section ^ to yield the the average probability for segmental motion 
as a function of a reduced temperature and lattice filling fraction for each of the architectures in Fig. ^ 



In Section 111 we combine our simulatiDu-cesults with the recently developed Born- Green- Yvon lattice model for the 
thermodynamic properties of polymeralj"li3 to investigate the mobility of the four polyolefins at given temperature 
and pressure. Our method of determining the friction coefflcients is not an absolute one, but using the probability for 
segmental motion of a linear chain (polyethylene in our case) at 413 K and atmospheric pressure as a reference value 
we are able to calculate relative values of the friction coefficients as a function of temperature and pressure. 

In order to compare our results with experimental data, we extract friction coefficients of the polyolefins of interest 



from experimental viscosity data as described in Section IV . For temperatures well above the glass transition temper- 



ature, we find that our approach gives a good qualitative representation of the variation of the friction coefficient with, 
temperature and chain architecture. We also compare the pressure dependence of the viscosity of polypropyleneEj 
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with our predictions and find very good agreement for low to moderate pressures. The results presented here as well 
as future directions of this work are discussed in Section M. 



II. SIMULATION OF LOCAL MOBILITY 

A central point of this work is the determination of the probability of segmental motion from a consideration of 
two short chain segments in a dense medium. To this end we perform an exact enumeration of all possible combined 
configurations and relative movements of two chain segments on a lattice. During the enumeration procedure we 
collect statistics on the characteristic parameters of each possible initial and final configuration and the connecting 
move. In a second step, these statistics are evaluated for conditions corresponding to different temperatures and 
densities. The advantage of this two-step procedure is that the time-consuming part, the exact enumerations, have 
to be performed only once to yield results that can be evaluated quickly for a variety of conditions. 

A. Procedure 

For each of two polymer molecules, referred to as chain one and chain two from now on, we consider a straight 
section composed of three repeat units with given (generally not identical) side group arrangements. The repeat unit 
in the middle is the section of interest in each case, while the attached units represent the rest of the (long) chains. As 
depicted in Fig. |l|, we employ repeat units with four carbon atoms in the backbone for all the polyolefins considered 
here. The simulation procedure is illustrated in Fig. |^. 

In a preparatory step, all possible conformations (side-group arrangements of the three-repeat-unit sections) of the 
chains are generated. The basic step in the enumeration, to be described below, is performed for given conformations 
and relative orientation of the chains. Keeping the conformation and orientation of chain one fixed, the basic step is 
repeated for all orientations of chain two. This sequence of steps is then repeated for different chain-two conformations 
until all conformations of chain two are exhausted. Then the whole sequence is repeated for all different chain-one 
conformations. This assures that we enumerate the results for all possible and distinct combined configurations of 
the chains. 

In a basic step, the section of interest of chain one is fixed to the origin of a simple cubic lattice and aligned with 
the z-axis. The total number nt and the coordinates of the non-bonded nearest neighbors (nn) sites are determined 
as is the maximum number Cm = 4st of possible contacts, where Sf denotes the number of lattice sites occupied by 
the section of interest. For a polyethylene (PE) chain, for example, the section of interest occupies Sf = 4 sites, has 
TT-t = 16 nearest neighbor sites and a maximum of Cm = 16 contacts, while for polyisobutylene (PIB) the values are 
Sf = 8, Tit — 24, and = 32. 

Next, an orientation for chain two is chosen, and the following procedure is repeated for each site of interest of 
chain two and each identified nearest neighbor site of chain one, where care has to be taken to avoid double counting. 
For the given orientation of chain two, contact is made between the chains by moving the currently considered site 
of chain two onto the currently considered nearest neighbor site of chain one. If the sections do not overlap, the 
combined configuration is accepted. It is evaluated by counting how many of the nt nearest neighbor sites of chain 
one are occupied by chain two (this number is called Oi) and by counting the number Ci of contacts between chains 
established in this way. The numbers (oi, Ci) characterize the static properties of the initial state. 

In order to determine the mobility of the segments, an attempt is made to displace chain one by one lattice site in 
each of the six directions, ±x, ±?/, and ±2:, in turn. If the attempt leads to overlap between the chains, it is counted 
as impossible. Otherwise, the number Sn of lattice sites newly occupied by the first section of interest is counted 
and the new combined configuration is characterized by determining the numbers o f and c/ of occupied nn sites and 
established contacts, respectively. The numbers {of,Cf) characterize the static properties of the final state, while the 
set (i, /, Sn) = {oi,Ci;of,Cf;Sn) characterizes the move. 

For the results presented here, the exact enumeration procedure described above was performed for PE, PEP, and 
PP. To avoid excessive computation times for PIB, we generated representative samples of one eighth of the single 
chain conformations and proceeded with those as described above. The different representative samples for the PIB 
conformations give essentially identical results for the properties presented in this work. The result of the simulations 
are the frequency of occurrence, n{i, /, Sn), of moves of type (i, /, Sn), as well as statistics on the type, (ofc, Ck), and 
frequency of occurrence, mk, of the combined configurations. 
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B. Evaluation 



In order to determine the probability P for segmental motion, we consider the probabilities P{i, f, Sn) for the 
different types of moves and then form the sum: 

P= J2 ^(*'/'^n). (2) 

The probability P(j, /, Sn) for a move of type (i, /, Sn) is expressed as 

P(^,/,.„) = !^Mi£:il^^, (3) 
2_^n{i,f, Sn) 

where Pi is the probability for the initial combined configuration to be of type {oi,Ci), Pae accounts for the energy 
difference AE — Ef — Ei between initial and final states, and is the probability that a sufficient number of 
contiguous sites is available to the moving segment in a lattice filled to a fraction (j). 

Both Pi and Pae involve the energy of a combined configuration of chain segments in a dense medium. As explained 
in detail in Refs. ^ and [T^, the energy E^. for a combined configuration characterized by (ofc,Cfc) is obtained from 

E, = eL+ ^'''^'"'^^'--''A , (4) 

V "t - Ofc / 

where e is the interaction energy between two molecular sites (e < 0) and where ^ = 20/(3 — (/)) is the contact density 
for infinitely long chains at a filling fraction of (j). The contributions to the energy Ek in Eq. (^) are due to contacts 
between chain one and chain two and between chain one and its randomly filled nearest neighbor sites, respectively. 
Please note that in the case of unbranched segments (PE), Eq. (^) implies Ek — eCmC for all k so that there are 
no energetically preferred configurations for PE. The unbranched segments thus serve as our reference system which 
allows the isolation of the effects of small scale architecture. 

The probability Pi for an initial combined configuration with energy Ei is proportional to the Boltzmann factor: 

P,^rr^e~^^^/Y,rr^ke-^''\ (5) 
fc 

where rui is the multiplicity of the combination (oi,Ci) and where /? ~ with temperature T and Boltzmann's 

constant k-Q. The effect of energetics on the probability of a move is described using the Metropolis form 

Pae-^^ ff as < ' 

where AS = Ef - Ei. 

The mobility of a chain segment is greatly reduced by the presence of the other chains in its surroundings. Consider, 
for the moment, monomers on a lattice and the attempt of a single particle to move from one site to a neighboring 
site. If all of its neighboring sites are occupied, the attempt will certainly fail. But even if a nearest neighbor site is 
available, the attempt may fail when another particle is headed for the same site. The only way to guarantee that 
an attempted move will be successful is to require that a neighboring site as well as its nearest neighbor sites (except 
for the one occupied by the particle under consideration) are empty. This is a total of six sites for monomers on a 
simple cubic lattice. Extending this reasoning to chain molecules, where each monomer is bonded (on average) to 
two monomers on neighboring sites, we require four empty sites in the neighborhood of each monomer involved in the 
move. In terms of our variables introduced above, a volume of 4sn is required in a move in which Sn sites are newly 
occupied. Assuming a random distribution of voids over the lattic&nthe probability of finding 4sn lattice sites among 
the Tit nearest neighbor sites of the segment of interest is given bjilS 

where (1 — (/)) is the fraction of empty sites. 

In Fig. 1^ we present simulation results for the probability of segmental motion, evaluated according to Eq. (^ with 
Eqs. (||) to (|^). The probabifity P is shown as a function of reduced temperature T/T* ~ k^Tje for a given filling 



3 



fraction 4> for the four architectures considered in this work. The effect of the increasing number of side groups on 
the mobihty is clearly visible. The linear chain (PE) has the highest probability of segmental motion followed by 
PEP, PP, and PIB which have one, two, and four side groups in the four-carbon backbone monomer, respectively. 
The insert shows the probability P as a function of reduced temperature for three different filling fractions for the PP 
architecture. As expected, the probability for segmental motion increases with reduced temperature and decreases 
with filling fraction. 



III. CALCULATION OF FRICTION COEFFICIENTS 



With a method to obtain the probability for segmental motion in hand, we are now in a position to address the 
monomelic friction coefficient C, oc P^^ . The proportionality constant between and C, is not easily determined. 
However, as pointed out earlier, our goal is to determine how the local architecture changes the mobility of branched 
chains compared to that of linear chains. Hence, we choose a reference state for the linear chain and express our 
results for the friction coefficients as the ratio 

where P^^i and Crof are the reference state values of the probability of segmental motion and the friction coefficient of 
the linear chain, respectively. The reference state can be chosen freely; in our case a temperature of Tiot=413.15 K 
and a pressure of pi.ef=0.1 MPa turn out to be convenient. In order to make contact with experimfiptal data we 
employ equations of state based on the recently developed Born- Green- Yvon (BGY) lattice modelliS'liil The BGY 
lattice model has three system-dependent parameters for a polymer melt, corresponding to the volume v per lattice 
site, the number r of sites occupied by each chain, and the interaction energy e between non-bonded nearest neighbors. 
For each of the polymers considered in this work, values for the system-dependent parameters have been determined 
from a comparison with experimental temperature-density-pressure datatJ and are summarized in Table ||. 

In Fig. ^ we present calculated values for the relative friction coefficient C/dcf as a function of temperature at a 
pressure of 0.1 MPa for the polyolefins considered in this work. Please note that the temperature variation here is 
much larger than that in Fig. ||. While Fig. || depicts the probability P oc C,~^ at constant 0, i.e. at constant density 
p = ip/rv, Fig. ^ shows constant pressure results. It is the temperature variation of the density at constant pressure 
that, through Eq. (Q), is responsible for the strong temperature dependence of C/Cef- 



IV. COMPARISON WITH EXPERIMENTAL DATA 



The friction coefficient ( is not a directly measured quantity but can be extracted from measurements of dynamic 
properties like the viscosity or the self-diffusion coefficient. The most direct access to C is through the Rouse viscosity 
77r described in the Introduction. With the aid of Eq. (|^) a friction coefficient per monomer can be defined as follows: 



^ _ a^rCr m [Na if 



(9) 



where N is the degree of polymerization and toq is the mass of a monomer which we take to be n x 14.03 g/mol, 
where n is the number of carbon atoms in the repeat unit depicted in Fig. |l|. In addition to the Rouse viscosities 
ryn/M, evaluation of Eq. requires values for Rq, the mean-squared end-to-end distance of the chains, and for p, 
the mass density of the melt. In this work, we employ experimental values for these properties at a temperature of 
413 K presented in a recent review by Fetters et aln and included in Table ||. 

The Rouse model describes directly the viscosity for melts of low molecular mass M <C Mo (for example, ~ 1000 
for PE and Mo « 7300 for PIBd). Unfortunately, measurements of the melt viscosities for short chains are not only 
scarce, but chain-end effects may have to be taken into account in theiii^valuation.D We therefore decided to turn to 
high molecular-weight viscosity data and an empirical scaling relationO to extract values-Jpr t/r/M. Motivated by 
the reptation modeO and by experience with experimental data, Graessley and Edwardalj suggested the following 
molecular mass dependence of viscosities in polymer melts; 



m 

M 



M 



M 
Mo 



(10) 
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with Mc = 2.2MoJli'0 In principle, Eq. (|o|) could be compared to experimental viscosity data at a given temperature 
to extract both rj^. and Mc- This, however, is not advisable since the results for rj^ and Mc will then strongly 
depend on the range of molecular weights for which viscosity data are available. Instead, we fix Mc = 2.2Mo, taking 
the experimental Me values at 413 K of Fetters et aln quoted in Table ||, and fit for rj^/M only. This procedure 
requires values for the_viacosity at 413 K, which we obtain by employing the temperature correlations provided with 
the experimental dataETQ to shift the viscosity values from the temperature of the measurements to 413 K. More 
information on the temperature correlations is provided below. In Fig. ^ we present the data (symbols) for the 
viscosity at 413 K obtained in this way as a function of molecular mass for the polyolefins considered in this work. 
The resulting correlations for the mass dependence of the viscosity are represented by the lines in Fig. ^ and are seen 
to give a satisfactory representation of the experimental data. The values for rj^^/M obtained in this way are included 
in Table ||. Inserting them into Eq. (^) and employing the values for toq, Rq, and p as discussed, we arrive at the 
values for the monomeric friction coefficient C at 413 K and atmospheric pressure presented in Tabic ||. The value of 
the friction coefficient for PE at 413 K and 0.1 MPa is the reference value for the experimental friction cocfhcients, 
and all further results will be presented as C/Cpe(413 K). 

In measurements of viscoelastic properties of polymers, it is customary to describe the temperature dependence of 
the viscosity by equations of the typed Vogel-Tammann-Fulcher (VTF) 



or Williams-Landel-Ferry (WLF) 



ln(r?(r)) = In(^) + (11) 
a{T - To) 



\og{v{T)) = log(7?(r,)) - '^f^^'^'l^ - (12) 



These equations are equivalent with constants related by C2 — Tg — Tq, C1C2 In(lO) — l/a, and ln(y4) ~ ln(?7(Ts)) — 
Ciln(lO). Here Tg is an arbitrary reference temperature while Tq indicates the temperature where the system is 
no longer able to relax to an equilibrium state in a finite amount of time. Since the dominant contribution to the 
temperature variation of the viscosity is due to the friction coefhcieiit, it is a reasonable approximation to assign 
the temperature dependence of the viscosity to the friction coefficient .111 Some of the experimental works quoted here 
provide slightly different temperature correlations for the viscosity of samples of different molar masses. When shifting 
the experimental viscosity data to the reference temperature of 413 K, we employed the correlations appropriate for 
the molar mass under consideration. For the following comparison with our work, on the other hand, we choose a 
representative correlation for each polyolefin and bring it into WLF form with a reference temperature of Ti.cf=413 K. 
The corresponding parameters Ci and C2 are included in Table ||. 

The temperature dependent friction coefhcients are now obtained from 

log(C(T)) = log(C(Tref)) - "^y^ (13) 

O2 + i — Jrcf 

with the C(Tlcf) values presented in Table |l|. In Fig. |^ we present the friction coefficients C{T) divided by the reference 
value Cpe(413 K) for the polyolefins considered in this work. The heavy lines in the graph indicate the temperature 
range in which experiments were performed and where the temperature correlations are expected to be most reliable. 
For each of the polyolefins, a strong increase in the friction coefficient is evident as the temperature is lowered. 
This increase is due to the slowing of the dynamics of the polymers as the glass transition is approached. The glass 
transition temperature Tg of polypropylene (aPP) is much higher (cf. Table O) than that of the other three polyolefins 
considered here. This is apparent in Fig. ^, where the friction coefficients of PE, PEP, and PIB have very similar 
temperature dependencies, while the aPP friction coefficient curve starts turning up at a much higher temperature 
and crosses the curve of the friction coefficient for PIB. 

In Fig. 1^ we present the friction coefficients extracted from experimental viscosity data together with those pre- 
dicted from our simulation procedure for temperatures well above the glass transition temperatures of the polymers. 
Comparing the predicted with the "experimental" curves we note that in both graphs the values of the friction co- 
efficients increase for a given temperature when going from PE, which has the lowest friction coefficient, over PEP, 
PP, to PIB. As noted earlier, this can be understood as a result of the different small scale architectures since the 
number of side groups in the repeat unit increases from PE (0) over PEP (1), PP (2) to PIB (4). The magnitude of 
the architecture effect is similar in the predicted and experimental friction coefficients. Furthermore, we note that the 
variation of the friction coefficients with temperature is of the same order of magnitude in the predicted and extracted 
curves. While the temperature range in Fig. |^ is well above the glass transition temperatures, lower temperatures 
are included in Figs. ^ and ^. As can be seen from these figures, the agreement between the simulation results and 
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the ^ values extracted from experimental data diminishes as the glass transition temperature is approached. This is 
because our simple approach focuses on individual segmental motion rather than cooperative effects, a point which 
will be discussed below. 

In Fig. H we compare the pressure dependence of the viscosity of polypropylene as predicted from our work with 
experimental datats along three isotherms. To separate temperature and pressure effects, we use the values of the 
viscosity at atmospheric pressure to scale the viscosities on each isotherm. For low to moderate pressures (<20 MPa) 
the agreement is excellent. As the pressure increases further, our work overestimates the viscosity. This isJikely due 
to our employing the BGY-lattice-model equation of state with parameters optimized for low pressures,E3 an issue 
which will be addressed in future work. 



V. DISCUSSION 

In this work, we presented an exact enumeration method for lattice simulations of chain segment mobility. The 
algorithm enumerates the attempted and successful moves for two short, straight sections of a polymer and is evaluated 
by taking relative frequency, energetics, and density effects into account. The result is the mean probability for 
segmental motion as a function of reduced temperature and the filling fraction of the lattice. We performed simulations 
for four different small-scale architectures obtaining results that show a sensible decrease in mobility with increasing 
density and number of side groups of a monomer. 

Combining these results with equations of state for the corresponding polyolefins, we deduce monomeric friction 
coefficients as a function of temperature and pressure. Our method is not an absolute one, but employing the friction 
coefficient of polyethylene at 413 K and atmospheric pressure as a reference value, we can predict the relative values 
of the friction coefficients as a function of temperature and pressure for the polyolefins considered in this work. If 
we are interested in the properties of a single polyolefin, we can employ the value of the viscosity for a particular 
temperature and pressure as a reference value and predict the relative variation of the viscosity with temperature 
and pressure from there. The same is true for other transport properties that depend in a simple way on the friction 
coefficient; this will allow us to investigate diffusion coefficients, for example, in future work. 

To compare our results with experimental data, we extracted Rouse viscosities and monomeric friction coefficients 
from high molecular mass viscosity data. Employing temperature correlations of the experimental data in the WLF 
form we obtain "experimental" monomeric friction coefficients at atmosperic pressure over a range of temperatures, 
which we scale by the value for PE at 413 K. The comparison of these extracted friction coefficients with the results 
from our new simulation method is encouraging: For temperatures well above the glass transition temperature the 
calculated probabilities give a good qualitative representation of the relative variation of the friction coefficient with 
temperature and monomer architecture. A comparison of calculated and experimental pressure variation of the 
viscosity of polypropylene along three isochores shows excellent results for lower pressures. 

In order to extend the range of validity of the present theory to temperatures closer to the glass transition, coop- 
erative effects in the dynamics will have to be taken into account in a more sophisticated way. In this first work with 
the new simulation method, we have assumed a random distribution of voids over the lattice and take the system 
to be in an equilibrium state before each attempted move. One way to improve on this approximation would be to 
employ an iterative approach in which the distribution of configurations after a round of attempted moves is used as 
the input distribution of configurations for the next round of moves. Finally, the simulation method introduced here 
was applied only to straight chain sections of polymers on a cubic lattice with a single site-site interaction strength e. 
It is, however, readily modified to include chain flexibility, chemical differences and realistic bond angles, which allows 
a large range of polymeric systems to be investigated in this way. We plan to extend the theory in these directions 
and are currently focusing on the effects of chain flexibility. 
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TABLE L BGY lattice-model parameters for the polyolefins considered in this work. Please note: The value r for the 
number of lattice sites per chain is proportional to the molar mass of the polymer. The values presented here correspond to a 
molar mass of AI = 170, 000 for each of the polyolefins. 



Polymer 


PE 


PEP 


aPP 


PIB 


e(J/mol) 


-1977.5 


-2000.0 


-2040.7 


-2208.1 


r 


14210.6 


14226.1 


14471.8 


13433.8 


i;(L/mol) 


0.013 


0.013 


0.013 


0.013 



TABLE II. Experimental parameters for the polyolefins considered in this work. The values correspond to atmospheric 
pressure and a temperature of 413.15 K. The glass transition temperature for polyethylene (PE) is an estimate based on results 
for ethylene-butene copolymersEj 



Polymer 


PE 


PEP 


aPP 


PIB 


R'o/M (A^mol/g) (Ref. |) 
p(g/cm«) (Ref.|) 
Me (Ref. 3) 


1.21 
0.785 

976 


0.834 
0.79 

2284 


0.67 
0.791 

4623 


0.57 
0.849 

7288 


viscosity data references 
temperature range T(K) 


m 

350-500 


i- 

248-443 


i- 

298-463 


1- 

298-473 


Inr-c'^R / Poise \ 
^ IT /g/mol^ 
log((;^/Poise cm) 


-4.1035 
-8.5554 


-3.6039 
-7.8001 


-3.3424 
-7.3648 


-2.7669 
-6.6388 


Ci 
Ca (K) 


2.018 
253 


3.565 
277 


3.101 
189 


4.684 
307 


T, (K) 
Tg references 


188 


211 
1 


268 



202 

SI 
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PE PEP PP 



PIB 



FIG. 1. United atom representation of the polyolefins considered in this work. Shown are the repeat units with four carbon 
atoms in the backbone for polyethylene (PE), an alternating copolymer of polyethylene and polypropylene (PEP), polypropylene 
(PP), and polyisobutylene (PIB). 
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FIG. 2. Illustration of the simulation procedure. The figure on the left shows an initial combined configuration of two PEP 
segments. There are Sf = 5 sites in the section of interest in each chain, indicated by the dark filled circles, and rit = 18 

identified nearest neighbor sites of chain one (indicated by open circles). The numbers o; and Ci indicate the number of nn sites 
occupied by chain two and the number of established contacts, respectively. The figure on the right shows the new combined 
configuration after chain one has been moved one lattice site to the right. The numbers Ok and Ck give occupied sites and 
established contacts, respectively. In the move, Sn = 5 sites were newly occupied by the first segment. 
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FIG. 3. Probability for segmental motion as a function of reduced temperature T /T* = ksT/e at constant filling fractions 
(j> determined from the simulation procedure and Eqs. (|2|) - (^. The graph shows results for the four different architectures 
(see Fig. ^) at a common filling fraction of = 0.8, while the insert displays results for three different filling fractions for the 
architecture of PP. 
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FIG. 4. Calculated friction coefficients C/Cref, as a function of temperature at constant pressure p = O.lMPa for the polyolefins 
considered in this work. The reference value is obtained for the linear chain (PE) at 413.15 K and 0.1 MPa. 
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FIG. 5. Viscosity of polyolefin melts as a function of molecular mass Afw The symbols represent experimental data adjusted 
to 413 K as explained in the text. The references for the data are Refs. for PE, | for PEP, for aPP, and | for PIB. The 

lines represent the empirical scaling law Eq. (|lo|), where Mc — 2.2Mc with Mc from Ref. ^ and where the values for rju/AI have 
been determined from a comparison with the experimental data. Values for these system-dependent parameters are included 
in Table 0. 
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FIG. 6. Temperature dependence of the friction coefficients ("(r)/("pE(413 K) extracted from experimental data. Tfie heavy 
hnes in the graph indicate the temperature range over which experiments were performed, cf. Table III 
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FIG. 8. Pressure dependence of the viscosity of polypropylene. For each of the three isotherms we present viscosities 
divided by their value at atmospheric pressure. |-3jhe lines represent results from our simulation procedure, the symbols indicate 
experimental data by Mattischek and SobczakJi3 
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